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NYO-61|86 
PliYSICS  AND 
MATHEMTICS 
PREFACE 

This  note  contains  an  expository  analysis  of  the  mathematical 
structure  of  the  theory  governing  interaction  of  magnetic  fields 
with  conducting  compressible  fluids.   One  of  its  major  aims  is  to 
emphasize  the  remarkably  close  pe4"u')llelism  between  this  theory 
and  ordinary  gas  dynamics.   It  is  shown  that  the  basic  equations 
governing  magneto-hydrodynamics  -  or  hydromagnetics  for  short  - 
have  essentially  the  same  mathematical  character  as  those  govern- 
ing gas  dynamics  and  that,  consequently,  essentially  the  same 
mathematical  methods  that  have  proved  successful  in  gas  dynamics 
can  be  employed.   This  fact  is  illustrated  by  a  detailed  descrip- 
tion of  the  hydromagnetic  analogues  of  shocks,  first  discovered  by 
Teller  and  de  Hoffman,  of  sound  waves,  including  the  Alfven  waves, 
and  of  simple  waves.   In  particular,  a  typical  example  is  presented 
which  serves  to  demonstrate  that  -  as  in  gas  dynamics  -  simple  one- 
dimensional  magneto-hydrodynamic  flow  problems  can  be  solved  with 
the  aid  of  shocks  and  simple  waves. 

The  work  presented  here  originated  in  connection  with  a 
Seminar  conducted  in  1951+  by  H.  Grad  at  New  York  University.   A 
preliminary  report  appeared  at  the  Los  Alamos  Scientific  Laboratory 
in  September  195I4-,  reissued  in  March  1957  as  Report  LAMS  2105o 

Some  of  the  problems  described  in  the  preliminary  report  were 
subsequently  treated  extensively  by  Bazer  and  by  Bazer  and  Ericson; 
their  results  are  referred  to  and  used  in  the  present  note.   Also 
included  is  an  appendix  by  K,  von  Hagenow,  in  which  certain  of  the 
results  on  simple  waves  are  shown  to  be  deduclble  from  the  trans- 
formation properties  of  the  differential  equations. 

The  relativistlc  analogue  of  the  problem  treated  here  was 
formulated  in  a  report  by  P.  Relchel  [12] j  see  also  Zuraino  [1?]. 
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NY0-6h.86 

NOTES  ON  MAf'NETO-HYDRODYNAMICS  -  iTrjT4P.Ea  VIII 
Nonlinear  Wave  Motion 

K,  0.  Priedrichs  and  H.  Kranzer 

(.      Introduction.   Physical  assumptions. 

The  equations  dealt  with  in  this  note  represent  a  part ic alar 
special  case  of  the  .full  set  of  raagneto-hydrodynamic  equations 
catalogued  in  MH-I,   This  specialization  is  achieved  through  a 
ntxmber  of  physical  assumptions  which  serve  to  reduce  the  mathe- 
matical complexities  inherent  in  these  equations  to  almost  manage- 
able proportions. 

Our  first  assumption  is  that  there  exists  a  scalar  fluid 
pressure  p  which  is  a  function  of  density  and  entropy.   This 
function  is  assumed  to  have  properties  usually  required  in  gas 
dynamics.   Moreover,  we  shall  ass\ame  in  our  discussion  that  heat 
conduction  and  viscosity  may  be  disregarded.   As  a  consequence 
various  types  of  gases  are  excluded  from  treatment,  such  as  geses 
in  which  the  mean  free  path  is  not  small  compared  with  the  signifi- 
cant dimensions  of  the  problem. 

We  also  shall  assume  that  the  flow  velocity  is  small  compared 
with  the  speed  of  light.   Pxirthermore,  we  assume  that  the  (mean) 
electric  charge  is  negligible,  so  that  the  medium  Is  essentially 
neutral,  and  that  the  displacement  current  may  be  neglected. 


*It  would  not  be  necessary  to  make  these  assiamptions,  A  strictly 
Lor entz- invariant  counterpart  of  the  treatment  given  in  this  report 
has  been  given  by  Reichel  [12]. 

-  h.  - 


On  the  other  hand,  we  do  not  assvune  that  the  flow  velocity 
is  small  compared  with  the  speed  of  sound.   In  other  words,  we 
assttme  the  fluid  to  be  compressible. 

Moreover,  we  shall  assume  the  electrical  conductivity  of  the 
fluid  to  be  infinite.   This  assumption  enables  us  to  express  the 
electric  field  in  terms  of  the  magnetic  field  and  the  flow  vector. 
Indeed,  from  equation  (26)  of  MS-IV  we  find 

(1)       E  =  B  X  u. 

It  is  not  our  intention  to  discuss  in  detail  the  significance 
of  these  various  physical  assumptions.   Instead,  we  want  to 
describe  some  of  the  mathematical  consequences  of  the  basic 
differential  equations  to  be  derived  from  these  assumptions, 

A  remark  might  be  made,  though,  about  the  assumption  of 
compressibility.   As  is  well  known,  compressibility  need  be 
taken  into  account  only  if  the  flow  speed  \u\    is  comparable  with 
the  sound  speed 


(2)       a  =  c^^^^  =  [dp/dp] 


1/2 


Moreover,  our  approach  will  in  general  be  worthwhile  only  if  the 
Alfven  speed 


(3)       °Alf  ~  '  ^^  *  ^^^   Chapter  II, 

is  comparable  with  the  sound  speed,  since  otherwise  hydromagne t i c 
and  compressibility  effects  could  be  separated.   If  c.-,-  is  inuch 
less  than  c     ,,  which  will  frequently  be  the  case,  it  might  be 
possible  to  separate  compressibility  effects  from  hydromagnetic 
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effects.   We  might  then,  for  example,  first  treat  compressibility 
effects  in  the  absence  of  hydromagnetic  effects,  and  then  consider 
hydromagnetic  effects  under  the  assumption  of  a  knovm  fluid 
density. 
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1,   Basic  Equations 

On  the  basis  of  all  the  assimptions  described  one  can  derive 
a  system  of  differential  equations  which  governs  the  flow  of  the 
fluid  and  the  changes  in  the  electromagnetic  field,  B  is  the 
magnetic  field  vector;  the  electric  field  vector  is  given  by  (1) 
in  terms  of  the  flow  velocity  Uc.  In  the  absence  of  displacement 
current,  the  current  per  unit  area,  J,  can  be  expressed  in  terms 
of  the  curl  of  3: 

(1.1)     ixJ  =   curl  B, 

The  pressure  p  =  p(p,S)  is  a  given  ftinction  of  density  p  and 
entropy  S, 

The  two  Maxwell  equations  which  do  not  involve  ciirreno  and 
charge  are  retained j  the;/  may  be  -written  ez 

A^     div  B   =  0 
o 


and 


B  +  curl  (B  X  u)  =  0. 


-"1 

The  second  term  in  A,  is  curl  E  by  (1), 

The  force  per  lonit  volume  which  enters  ."ewton's  second  lav: 
consists  of  the  "Lorentz  force"   -B  x  J  and  the  pressure  gradient 
"VP*   '•■sing  (1.1),  we  may  write  this  law  in  the  form 

Ag     pu  +  p(u«  y)u  +Vp  ■*■  ^i-'-B  X  curl  B  =  0. 

The  continuity  equation  of  fluid  dynamics  is  retained; 

A  2  p  +  c  i  V  ( pu )   =  0 . 

Finally,  vie    add  the  law  that  the  entropy  per  unit  riass,,  S,  is 
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c&rried  unchanged  by  the  particles: 
A,      S  +  (u.V)S  =   0, 

These  equations,  essentially  formulated  by  Lundqulst,  1952, 
[3],  and  sometimes  referred  to  as  "Lundquist  equations",  will  be 
the  basis  of  our  discussion.  They  include  as  a  special  case  the 
two-dimensional  problems  treated  in  MH-VII. 

It  ivould  be  sufficient  to  require  equation  A   to  hold  only 
at  an  initial  time;  as  is  well  known,  it  then  follovjs  from  A^  that 
this  equation  is  satisfied  at  all  times. 

Equations  A-.  to  Aj  are  a  system  of  eight  non-linear  partial 
differential  eauations  of  first  order  for  the  eight  quantities 
B  ,  B  ,  B  ,  u  ,  u  ,  u  ,  p,  and  S.   The  first  fact  we  emphasize  is 
that  these  equations  are  hyperbolic.   Specifically,  they  belong 
to  the  special  class  of  symmetric  hyperbolic  equations,  which  is 
particularly  well  understood  mathematically;  see  [II4.].   These 
equations  share  this  property  with  Maxwell's  equations,  with  the 
equations  of  elasticity,  and  with  the  equations  of  gas  dynamics. 
Moreover,  they  share  with  the  gas-dynamical  equations  the  further 
property  of  being  nonlinear  with  coefficients  which  involve  the 
dependent  but  not  the  independent  variables.   It  is  because  of  the 
latter  property,  in  addition  to  the  hyperbolic  character,  that  the 
same  methods  can  be  applied  to  the  equations  A,  to  Ai  that  have 
proved  successful  in  gas  dynamics. 

The  fact  that  these  equations  of  magneto-hydrodynamics  belong 
to  a  class  of  equations  which  are  most  consistent  in  their  mathe- 
matical character  supports  the  confidence  in  the  physical  consis- 
tency of  the  various  assumptions  from  which  they  were  derived, 
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As  is  well  known,  physical  processes  governed  by  hyperbolic 
equations  have  the  property  that  distirrbances  sve   propagated 
with  finite  speed.   Thus,  in  ordinary  compressible  fluids, 
disturbances  travel  with  the  speed  of  sound  relative  to  the 
motion  of  the  fluid.   Hence  magneto-hydrodynamic  disturbances 
also  travel  with  finite  speed.   However,  in  contrast  to  gas 
dynamics,  there  are  three  "sound  speeds,"   Accordingly,  there 
are  three  modes  of  propagation  in  each  direction.   Moreover, 
these  speeds  depend  on  the  directionj  specif ically^,  on  the  angle 
between  the  direction  of  propagation  and  the  direction  of  the 
magnetic  field.   These  remarkable  facts  were  first  discovered 
by  Herlofson  [3]  and  by  van  de  Hulst  [I|-],  in  connection  with 
sinusoidal  linearized  wave  motion.   Our  first  task  will  be  to 
find  these  speeds. 
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2o   Characteristic  Manifolds  and  Propagation  of  Dlatxirbances 

Characteristic  manifolds  -  three  dimensional  manifolds  in 
(x,y,z,t )-space  -  associated  with  a  differential  equation  may  be 
defined  in  many  different  waysj  cf„  [li^^l^lo   Instead  of  giving 
a  precise  definition^  it  is  sufficient  in  the  present  context  to 
say  that  solutions  of  the  differential  equation  may  possess  "small" 
discontinuities  only  on  certain  manifolds,  and  that  such  manifolds 
are  called  characteristic o   We  may  consider  such  a  manifold  as 
being  swept  out  by  siorfaces  sS  =fS{t)    in  (x,y3z)-space5  the 
motion  of  these  surfaces  will  then  also  be  called  "characteristic" o 
In  a  process  described  by  a  solution  of  the  differential  equation, 
therefore,  a  small  discontinuity  or  "disturbance"  present  on  a 
surface  JS*  ( t  )  at  an  initial  time  t  may  at  later  times  be  present 
only  on  surfaces  i^i  {t)   which  move  characteristicallyo   Such  a 
moving  discontinuity  will  be  called  a  "distiirbance  wave"  or  simply 


a  "wave",, 


We  introduce  the  normal  vector  n  of  xinit  length  at  each  point 
of  the  surface  r^J  ( t )  and  characterize  the  motion  of  the  surface 
(2^(t)  by  its  velocity  c  .  in  the  normal  direction  at  each  of  its 
points o   It  is  convenient  to  Introduce  the  normal  component 

{2ol)      u  -   nou 

©f  the  flow  velocity  at  this  point  and  to  write  the  characteristic 
velocity  c  ,  there  as  the  sum 

^2''2)     c^^     -  u^t   C5   • 

by  convention  we  always  choose  c  >  0, 
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Thus  Ic  is  the  normal  component  of  the  characteristic  velocity- 
relative  to  the  flow  velocity. 

In  order  to  find  the  possible  values  of  c  one  may  first  set 
up  the  relations  between  the  possible  discontinuities  5B,  5Uj  5p, 
6S  of  the  quantities  B,  u,  p,  S  on  the  surface  \!^  ( t ) .   Using  the 
formalism  of  the  theory  of  characteristics,  cf,  [li|-],  the  following 
relations  are  foijnd; 

B,         Tc5B  +  B5u   -  B  5u  =  0, 
1        +         n    n       ' 


Bg        +pc6u  +  a  n5p  +  p."-'-n(B<'6B)  -  |j,'""^B  5B  =  0, 


B 


3        +c5p   +   p5u^  =  0, 


B,.        +c5S  =  0 


Here  a  is  the  speed  of  sound,  given  by  (2);  furthermore 


(2o3)     5u  =  n»6u. 


and 

i2ok)  B^  =  n»B 


is  the  normal  component  of  the  magnetic  field. 

The  determinant  of  this  system  of  eight  homogeneous  equations 
for  the  eight  quantities  53,  5u,  5p,  6S  is  found  to  be 

(2.5)      det  (B)  =  pc^(pc^-  B^/[i)lpc^-    ( pa^+B'-/ti)c^+  ^^^n^^^  }  " 

The  characteristic  velocities  +c  are  obviously  the  roots  of  the 

equation 

(2o6)     det  (B)  =  0, 

'  "   11    " 


All  eight  roots  of  this  equation  are  seen  to  be  real^  in 
accordance  with  the  fact  that  the  differential  equations  are 
hyperbolic o 
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3 o   Fast  and  Slow  Disturbance  Waves 

We  shall  first  discuss  the  roots  of  the  last  factor  of 
det  (B)o   The  condition  that  this  factor  vanishes  can  be  written 
in  the  form 

(301)  c^Cpc^  -  B^ii)   ^  a^Cpc^  =  B^/|i.) 

or  in  the  form 

(302)  (c^  -  a2)(pc2  "   B^/ji)  -  c^CbV^J^  "  B^/V), 

The  larger  and  the  smaller  of  the  roots  c  >  0  of  this  equation 
will  be  denoted  respectively  by  c„   .  and  c  .   » 

Prom  equation  (3o2)  one  immediately  deduces  the  inequalities 


(^°3)  ^slow^  ^^  ^fast 

and 

(3oh.)  c  T,   <  b  <  c^   .  , 

^  ^  slow  °  n  "  fast 


Here 


(3o5)     b^  -  [Byii,p] 


1/2 


is  the  Alfven  velocity^  cfo  iDp   with  the  magnetic  field  vector 
B  replaced  by  its  normal  component  B  o   The  sovaid  speed  a  is  as 
given  by  (2)o   An  equality  sign  can  hold  in  relations  (3o3)  and 
(3oli)  only  if  B  =  B  n^  so  that  the  right  member  of  (3o2)  vanishes, 
In  this  case  one  of  the  two  speeds  equals  Ep  the  othiar  equals  b  » 

The  possible  disturbances^  ioeos  the  solutions  of  equations 
B,  to  B.  associated  with  c  =  ^gi^M  or  o  -  ^^fa^t'  ^^®  found  to  be 
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(3o6)     5B  =  kpc^(B  -  B^n), 


5u  -  ^kc(!a"^B^^B  °  pc^n) 


5p  =  kp(pc^  =.  B^/^), 


53  =  0. 


Here  k  Is  any  ntxmber  ^  Oo 

It  is  to  be  noted  that  5B  has  a  tangential  direction  so  that 

(307)  5B^  -  n^SB  -   0„ 

The  disturbance  5p  vanishes  when  c  =  b,  <,   The  distxarbance  5u  may- 
be written  in  the  form 

(308)  5u  =   f.Ciipcl^'^B^SB  +  p°^c(5p)n, 

which   is   on  occasion  useful o      We   also  note   dovm  the  relation 

(309)  5(p  +  B^a^t)   s  a^5p  +  |i°^°5B  ^  kpc^(pc^=  B^/ti,}, 

which  follows  from  {3o6)  and  will  prove  useful  later  oHo 

Particular  attention  should  be  paid  to  the  cases  where  the 
normal  n  is  parallel  or  perpendicular  to  the  magnetic  field  Bo 

In  the  first  casep  B  ■=  B  n,  one  of  the  roots  c  agrees  with  the 
sound  speed  ao   Unless  a  =  b  =  bj,  formulas  (3o6)  remain  valid 
for  this  rooto   They  must  in  any  case  be  modified  for  the  other 
root  c  =  bg  since  the  latter  root  agrees  with  a  root  of  another 
factor  of  det  {B)o 

In  the  second  casCp  B  =  0^,  the  fast  speed  is  given  by  the 

=  lij.  = 


noteworthy  formula 

Formulas  (3<,6)  remain  valid  in  this  caseo   The  other  root 
(3.11)    c^^^^  =  0, 

however,  agrees  with  the  root  of  another  factor  of  det  (B). 
Hence  for  this  root  formulas  (3o6)  must  again  be  modified. 

The  needed  modification  of  formulas  (3,6)  in  these  two 
cases  will  be  described  later  (^cf,  (i;o6)  and  {5.l4-))o 

In  order  to  illustrate  disturbances  of  the  type  treated  in 
this  section,  consider  segments  of  plane  wave  fronts,  fast  and 
slow,  which  travel  in  the  direction  of  the  normal  n  after  having 
passed  through  the  origin  at  the  time  t  =  Oo   At  a  time  t  >  0 
these  fronts  pass  through  the  points  c  t  n  with  c  =  c„   .  and 
c  =  c  ,   o   The  locus  of  these  points  is  shovm  in  Figure  1. 
Also  of  significance  is  the  envelope  of  these  fronts,  which  is 
shown  in  Figiire  lAo   This  envelope  gives  the  position  at  the  time 
t  >  0  of  the  expanding  wave  front  resulting  from  a  point  distur- 
bance at  the  origin  at  t  =  0 „   (That  is,  it  represents  the 
limiting  form  for  large  t  of  the  wave  fronts  resulting  from  an 
initial  disturbance  in  a  res'ion  of  finite  extent — for  small  t 
such  wave  fronts  may  look  quite  differento)   In  these  figures  we 
have  assumed  c.,^  <  c     ,,  so  that  the  slow  speed  c  ,    agrees 
with  the  Alfv^n  speed  for  B  |  |  n.   This  seems  to  be  the  more 
frequent  situationo   We  also  have  drawn  wave  fronts  of  a  third, 
intermediate,  type  which  we  are  going  to  describe  now, 
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lj.o   Transverse  Waves  and  Alfven  Waves 

An  intermediate  wave  belongs  to  the  root 

of  equation  (2o5)o   As  seen  from  Oolj.),  the  speed  b   lies  between 
c  ,    and  c~   .  except  if  B  i|n  or  B  J_  no   The  possible  discon- 
tinuities associated  with  this  wave,  given  as  the  solutions  of 
equations  (B),  are 
ihoZ)  5B  =  +  kpcB  X  n, 

6u  =  ku""'"B  B  X  n„ 

— r-     jj        J- 

5p  -  5S  -  Oj, 

with  an  arbitrary  constant  k.   In  this  case,  therefore ,  the  dis= 
ttirbances  are  tangential  to  the  wave  front  and  perpendicular  to 
the  magnetic  fieldo   The  relation 

(i|o3)      6B^   ^  SB-'SB   ^  0, 

Implied  by  (l|.o2)j,  shows  that  the  magnetic  field  xindergoes  a 
rotationo   This  intermediate  type  of  disturbance  wave  will  also 
be  referred  to  as  a  "transverse  wave" « 

At  this  place  we  may  interpose  a  remark  about  the  Alfven  wave 

in  an  incompressible  fluid o   The  conditions  on  the  possible  disturbances 
in  this  case  are  obtained  from  conditions  (B)  by  setting  a  5p  =  5p 
in  Bp  and  6p  -  0  in  B^o   In  addition  to  the  double  root  c  =  0,  one 
finds  that  c  -  b   is  a  double  root.   The  corresponding  distvirbances 
are 
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ikoh)  5B  ^  kB^n* 


6u  '-   Xkb  n* 

^   n 


f>p  =  -k|i"-^B^(B°n'""} 


5S  -  Op 


■where  n'"'  is  an  arbitrary  unit  vector  perpendicular  to  n,  and  k 
is  an  arbitrary  numbero   Evidently,  these  possible  disturbances 
form  a  two-dimensional  nianifold,  in  accordance  with  the  fact  that 
c  =  b   is  a  double  rooto   Relations  (L|.oi4.)  imply  the  important 
relation 

ikoS)  fi(p  +  bV?|a)   -  QP  +  p."^B^^6B  ==  0, 

which  expresses  the  fact  that  the  sum  of  fluid  presstire  p  and 

"magnetic  pressure"  B''/2|.t  is  continuous  across  the  surface  !^(t)o 
(The  notion  of  magnetic  pressure  will  be  discussed  in  Section  60) 

One  might  say  that  the  Alfven  wave"  results  if  the  sound  speed  a 

and  hence  the  fast  distiii^l")ance  speed  c-.p      .    become  infinite  while 

the  slow  vjave  speed  c  ..,         coalesces  with  the  Alfven  wave  speed 

"■  s  low 

For  a  compressible  fluid,  we  have  already  mentioned  an 
analogous  case  in  which  the  slow  (or  possibly  the  fast)  disttir- 
bance  speed  coalesces  with  che  Alfven  wave  speedy  roameljo  the  case 
of  a  disturbance  wave  traveling  parallel  to  the  magnetic  fields 

B   "  B  Ho 
n 


We  use  this  terra  for  distujr'bance  waves  in  any  direction  in  an 
incompressible  medium  although  Alfven  described  onljf  waves  travel- 
ing in  the  direction  of  B<, 
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2    2     2 
Again  c  =  b   =  b   is  a  double  root  of  equation  (2o6)  and  there 

again  exists  a  two -parameter  family  of  disturbances,  namely 
(ij.o6)     5B  =  kB^n-"-, 

5u  =  ^  kbn'"', 

6p  -  5S  -  Oo 

Here  n'""  and  k  have  the  same  meanings  as  in  equation  (l4.oi4-)o   Once 
more  we  have 

(l|o7)      5(p  +  B^S^i)  ^  ii°°-'-B°6B  =   0„ 

The  disturbances  (I;o6)  represent,  in  a  real  sense,,  the 
transition  between  a  slow  (or  fast)  wave  and  a  transverse  wave,, 
Indeedj  it  is  easy  to  see  that  they  are  obtained  as  the  limit  of 

equations  (l+<,2)  as  B  approaches  the  normal  directlong  and  also  as 

2     2     2 
the  limit  of  equations  (3o6)  as  B  <=>  B  n  and  c   =■>  b   =>  b  « 

The  fast,  slow,  and  Intermediate  waves  described  in  Sections  3 

and  l\.   represent  the  three  modes  of  wave  propagation  referred  to  at 

the  end  of  Section  lo   Since  each  type  of  wave  may  move  in  the 

direction  of  n  or  of  -n,  they  correspond  to  six  roots  Co   The 

propagation  of  disturbances  associated  with  the  remaining  roots 

c  ==  0  will  not  be  referred  to  as  wave  raotiono 


18 


5o   Contact  Distixr-bances 

About  the  remaining  root   c  =  0  only  a  few  remarks  need  be 
madeo   This  root  has  the  multiplicity  2  unless  B   =  0,  in  which 
case  the  multiplicity  is  6o   We  confine  ourselves  for  the  present 
to  the  case  B  f^  0.   The  relations  (B)  then  possess  only  the 
solutions 
(5ol)     5B  =  kn 

5u  -■      Og  5p  =  Oj 

5S  =  k^, 

with   aj?bitrary  k   and  k..  ,      We  may   just   as  well    set  k  =  0,    or, 
equivalentljj 

(5c2)  5B^     =      0„ 

It  is  consistent  to  do  so,  for,  as  could  be  shown,  this  condition 
is  satisfied  on  every  surface  j^(t)  if  it  is  satisfied  initially 

on  ^  ( t,J  o 
o 

The  remaining  possibility  of  an  entropy  disturbance  corres- 
ponds to  a  contact  discontinuityo   Since  the  discontinuity  is 
small,  we  prefer  to  use  the  term  "contact  disturbance". 

In  striking  contrast  to  the  situation  in  gas  dynamics  a  contact 
discontinuity  in  a  hydromapnetic  fluid  does  not  permit  a  discon- 
tinuity in  the  tangential  component  of  the  velocity  (provided  B  j^  0) 
This  remarkable  fact  will  play  a  considerable  role  in  the  discus- 
sions of  wave  motions  given  in  Section  13« 

Suppose  now  that  B  i£  perpendicular  to  n: 

B^  =   0, 
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Then  we  have 

so  that  c  =  0  is  a  sextuple  root  of  (2o6,)o   The  corresponding 
six-parameter  family  of  possible  disturbances  may  be  written 
as  follows g 

(5ol|)      SB  =   K, 

5u  =   kn's 

^1  =2 
5p  ~  =.^i  a  K'B^ 

5S  -   k,  o 

Here  K  is  an  arbitrary  vector ^  n*""  an  arbitrary  unit  vector 
perpendicular  to  n^,  and  k  and  k=,  arbitrary  numbers o   In  other 
wordSj,  all  the  disturbances  are  arbitrary,  except  that  they  must 
«^  satisfy 

(5o5)     6(p  +  B^/2p.)   -  Oo 

Equations  (5oI|)  do  not  form  as  clear-cut  a  "transition  case" 
as  equations  {\ob)  r.    since  performing  the  appropriate  limiting 
processes  on  equations  (3o6)p  (l|o2)j,  or  (5ol)  will  lead  only  to 
special  cases  of  (5ol4.)o 


20 


6o   Conservation  Laws 

In  this  section  we  shall  discuss  the  conservation  form  of 
the  Lundqulst  eauatlons,  partly  as  preparation  for  the  discussion 
of  shocks o 

It  is  customary  to  say  that  a  system  of  partial  differential 
equations  has  "conservation  form"  if  each  equation  consists  of 
the  sura  of  derivatives  of  functions  of  the  unknown  quantities. 
The  reason  for  this  expression  is  that  the  laws  of  conservation 
of  masS;,  momentum  and  energy  are  of  this  formo   On  the  other  hand 
the  possibility  of  writing  the  equations  in  this  form  is  the  condi- 
tion for  the  possibility  of  setting  up  shock  relations. 

Equation  A   (see  Section  1)  evidently  has  conservation  form, 

Cq         V»3   =   Oo 

This  same  is  true  of  equation  A,,  but  we  prefer  to  write  it  in 
the  form 

C^        B  +  V'uB  "      V-Bu  =  0. 

Here,  and  correspondingly  in  the  following,  the  inner  product  in 
V  'UB  involves  only  V  and  u  while  the  differentiation  ^  applies 
to  the  product  xjB.   Thus,  the  x-coraponent  of  \7  "UB  is 

(  V"UB)    -  #-  {u  B  )  +  J-  (u  B  )  +  -I-  (u  B  ). 
^  X    6x   XX    oy   y  X    oz   z  x 

Similarly,  the  x~component  of  V'Bu  is 

(  V°Bu)   =  ^  (B  u  )  +  J-  (B  u  )  +  I-  (B  u  ). 
^     X    dx   X  X    oy   y  x    oz   z  x 
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Equation  Ap  can  be  written  in  the  form 

Cg        ^  (pu)  +  V"u(pu)  +  7p  +  VCbVsm-)  "  V°Bix"^B  =  0, 

which  expresses  the  law  of  conservation  of  momentxario   Suppose  the 
magnetic  field  B  is  such  that  the  term  V-B^"  B  V8nishes|  this  will 
be  the  case  under  various  symmetry  conditions,  cfo  [9K   Then  the 

conservation  law  Cp  has  essentially  the  same  form  as  in  gas 

p 

dynamics  except  that  the  expression  p  +  B  /Z^i.  takes  the  place  of 

the  pressure,.   It  is  primarily  because  of  this  fact  that  the  term 
B  /2jx  Is  referred  to  as  "magnetic  pressure",, 

The  law  of  conservation  of  mass^  A-.,-,  is  given  In  conservation 
form 

C,        p  +   Vm'pu)  ■=  Oo 

Equation  A.  ,,  however^  which  describes  the  transport  of  entropy^ 
Is  to  be  replaced  by  the  law  of  conservation  of  energy,  which. 
assumes  the  form 

\  ^  ^\^^   "^  P^  ^^   BV2|iy+  V°u(|pu2  +  pe  +  p  +  B^/(x) 

Here  e  is  the  internal  energy  per  unit  mass  of  the  fluids  which 
may  be  considered  as  a  function  of  density  and  entropy  and  is 
characterized  by  the  relation 

(6„1)     de  ^  TdS   =   pd(p°-^)^ 

in  which  T  is  the  temperature.   It  is  to  be  noted  that  the 
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expression  B  /2p,  in  the  first  term  of  C,  plays  the  role  of 
"magnetic  energy  per  unit  voliime" ,  while  the  expression  B  /^ 
in  the  second  term  of  Ci  plays  the  role  of  "magnetic  enthalpy  per 
unit  volume" o   It  should  also  be  noted  that  the  expression 

ubVm-   -  ^"■''B(B»u) 

occurring  in  C.  may  be  written  in  the  form 

{6o2)      ti"-^[viB^  -  B(B«u)]  =  ^"■'-[E  X  B]   =  E  X  H 

by  virtue  of  formula  (1).   This  term  is  thus  recognized  as 
Poynting's  energy  flux  per  unit  area. 
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'7o   Shocks 

As  is  well  knowrij  continuous  gas  dynamical  motions  will 
break  down  at  some  time  if  they  Involve  a  compressiono   The  same 
must  be  expected  to  happen  in  hydro-magnetic  motions  if  the 
compressibility  of  the  fluid  cannot  be  neglectedo   Mathematically 
speaking,  the  solution  of  the  differential  equations  ceases  to 
exist  beyond  the  time  of  breakdowno   Physically  speaking^  the 
phenomenon  is  no  longer  governed  by  the  differential  equations 
from  that  time  ono   Actually j,  discontinuities^  ioeoj,  shocks,  will 
appear  o 

As  in  gas  dynamics ,  one  assumes  that  the  quantities  on  both 
sides  of  a  shock  front  are  governed  by  the  laws  of  conservationo 
The  shock  conservation  laws  can  then  be  derived  from  the  conserva- , 
tion  form  (C)  of  the  differential  equations  by  a  formalism  which 
is  quite  analogous  to  the  formalism  by  which  the  equations  (B)  for 
the  disturbances  are  derived  from  the  original  differential 
equations  (A)<,   We  denote  by  n  the  normal  vector  at  any  point  of 
the  shock  front  and  by 

[Q]   =  Q.  -  Q„ 

1    o 

the  jump  of  any  quantity  Q  across  the  shock  frontg  here  Q.^  is  the 

value  on  the  side  toward  which  the  normal  n  points  and  Q   is  the 

o 

value  on  the  other  side„   Purtherp  we  denote  by  U  the  velocity  of 
the  shock  front  in  the  normal  directiono   The  recipe  then  requires 
one  to  replace  the  symbol  ■yr-  in  equations  (C)  by  •=U[^.ol  and  the 
symbol  V  by  n[  <> » ,  ]  o   The  result  is  the  following  set  of  equationso 
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^1  t(u^   -   U)B   -   B^u]    =   0, 

°2  f(^n   -   U^P^  +   <P   +  BV2n)n   -   ii'\B]    =  0, 

D3  [(u^   -   U)p]    =   0,  i 

These  are  the  relations  which  connect  the  values  of  the  cu.nti- 
ties  B,  u.  p.  S  on  one  side  of  a  shock  front,  or,  more  generally, 
of  a  discontinuity  surface,  with  the  values  of  these  ouantities 
on  the  other  side  and  with  the  speed  U  of  the  sur-face. 

These  relations  were  derived  by  de  Hoffmann  and  Teller  [2] 
m  1950„   These  authors  set  up  the  conservation  laws  for  shocks 
directly  without  relating  them  to  differential  eouatlons.   Also, 
they  derived  the  equations  in  a  Lorentz  invariant  form  and  only 
afterwards  derived  equations  (D)  as  the  non-relativistic  approxi- 
mation.  The  direct  non-relativistic  derivation  was  given  by  Ltst    ■ 
[Tlo 

The  jump  conditions  (D)  are  frequently  supplemented  by  the 
statement  that  there  is  a  "sheet  current"  flowing  along  the  dis- 
continuity surface  and  that  the  value  r   of  this  current  per  unit 
width  is  given  by  the  relation 

(*)       \if"   =  n  X  [B] 

in  accordance  with  relation  (l.l).   Although  this  statement  is  of 
great  significance  for  the  description  of  the  physical  phenomena 
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involved,  it  need  not  be  taken  into  account  in  the  mathematical 
analysis  of  the  discontinuity  conditiono 

We  will  speak  of  a  shock  if  fluid  crosses  the  front o  "u  -  V  ^   Oo 
In  that  case  vie  assume  the  direction  of  the  normal  vector  n  so 
chosen  that  fluid  crosses  in  the  direction  of  the  normals 

(701)  u^  =  U  >  Oo 

If  u  -=  XJ  —  0,  so  that  no  fluid  crosses  ^  we  speak  of  a  contact 
discontinuitjo 

We  shall  present  an  analysis  of  the  possible  types  of  shocks 
and  contact  discontinuities  which  is  completely  analogous  to  the 
analysis  of  disturbance  waves  given  in  the  preceding  sections o 
Before  doing  this,  however j,  we  make  an  important  observation  of 
a  general  natioreo 

Gas  dynamical  shocks  involve  a  rise  in  pressure  and  densityj 
since  the  entropy  increases  across  a  shocko   For  hydromagnetic 
shocks  we  may  state  similarlys   If  the  entropy  Increases  across  a 
shock  fronts,  pressxire  and  density  also  increase o   (The  qualifica- 
tion is  necessary  in  this  case^,  for,  as  we  shall  see^,  there  are 
shocks  which  do  not  involve  changes  in  entropy^  density j,  and 
pressure  at  alio )   The  proof  of  the  statement  follows  from  the 
identity 

(702)  [e  +  pP°^"]  -  =^[p"^][B]2/l4A, 

in  which  p  is  the  mean  value 

(703)  p  -  ^(Pq  +  p-l)o 
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This  remarkable  identity,  first  given  by  Liist,  could  be  derived 
by  forming  a  linear  combination  of  relations  (D)  which  is  similar 
to  the  linear*  combination  which  expresses  the  entropy  relation 
Aj  in  terms  of  the  relations  (C),  cf.  Section  6.   The  left  member 
of  (7o2)  is  the  "Hugoniot  function",  which  vanishes  for  gas 

dynamical  shocks.   The  right  member  involves  the  drop  -[p   ]  in 

-1 
specific  volume  p    and  the  square  of  the  Jump  of  the  magnetic 

fieldo   From  the  known  properties  of  the  Hugoniot  function 

—1 
(cfo  [12]),  one  can  show  that  it  has  the  same  sign  as  -[p   ] 

only  if  it  has  the  same  sign  as  [S],   The  statement  made  above 

then  follows. 

Prom  relation  (7o2)  one  can  also  derive  the  fact  that  the 
increase  of  entropy  across  a  shock  is  of  the  third  order  in  p 
and  Bo 

In  order  to  establish  the  analogy  between  the  shock  relations 
(D)  and  relations  (B)  characterizing  disturbances  it  is  convenient 
to  introduce  the  notion  of  mean  value 

Q  =  |  (Q^  +  Q-L^ 

of  any  ouantity  Q  and  to  use  the  formula 

[PQ]  =  P[Q]  +  [P]Q. 

It  is  also  convenient  to  introduce  the  specific  volume 

instead  of  p,  and  the  flux 
(7o5)     m  =-  p(u^^  -  U) 
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instead  of  u  «   Note  that^  by  D^,  the  flux  is  the  same  on  both 
sides  of  the  shock  fronto   AlsOji  D  permits  us  to  take  B^  as  a 


n 


constanto   Relations  D-.  to  D,  can  now  be  written  as 

E^        m?'[B]  +  B[u^]  "  B^[ul  =  0, 

E2        m[u]  +  [pin  4  n°^nB"[B]  =  ^i°^B^[B]  =  0, 

E^        m[r]    -    [u^J  "  0„ 

These  equations  evidently  correspond  precisely  to  equations  B, 
to  B^  if  one  lets 

m,  f  J,  ■=['f]   tp],B  correspond  to  ^pc,  p"  ,  p  a  ,  B 

and  lets 

in] g    [V] g    [B]  correspond  to  5u,  -p   Sp,  6B0 

The  analogue  of  relation  B.  is  relation  (7o2)j  it  may  be  disre= 
garded  in  the  present  context o 

Prom  this  analogy^^  or  by  direct  computations  one  finds  that 
the  determinant  of  the  system  E^  to  E-j  is 


^sj'' 


'^m^T^n,^^     U^  /,.   \  f   r  mh    +      (  T     m=lrr^1„     U^  /,   v      2 


(7o6)  det    IE)    -     r'^mCrm"^-  By|i)(rm'^  +    (TT  [tl'-'ipl^   BVli) 


n" 


m 


•ixrhv}  ^Ih) 


The   equation 

(7e7i  det    (E)      ^     0 

is  an  equation  for  the  flvix  m,  but  it  may  just  as  well  be  consid= 
ered  an  equation  for  the  shock  velocity 

(7o8)     u  ==  u^  =.  mr, 
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cf.  (7.5). 

It  is  clear  from  (7.6)  that  there  are  fast,  slow,  and 
Intermediate  shocks  and  that  the  relationship  between  them  is 
the  same  as  that  between  the  corresponding  distiorbance  waves. 


-  29  - 


8o      Fast   and   Slow   Shocks 

The   fliix  ra  of   a  fast   or   a   slow   shock  satisfies   the   equation 

(8ol)  m^Crm^   -  B^/y.)   =  =.[r]°^[p]  (  rm^=  b|/^), 

which  expresses  the  condition  that  the  last  factor  of  det  (E) 
vanlshesi  cfo  (3ol)o   In  analogy  with  (3o3)  and  (3ol+)  we  have 

and 

<8o3)     ^slow^^/^^^  ^4ast° 

In  analogy  with  (306)5,  the  relations  between  the  possible  jtmips 
across  the  shock  front  are 

(8oi^.)      [B]  -  k^m^CB  =  B^n), 
(8o5)     [u]  ^  k^^mC^i^^B^B  »rm^n), 
(8„6)     [r]  -  =k^(Tin^  =  ^n^^^» 

with  an  arbitrary  constant  k,  (instead  of  Tk)o 

The  cases  of  a  "parallel"  shocks  B  -  B  n^,  and  of  a  "perpen= 
dicular"  shocks  B  ~  0^  were  discussed  in  great  detail  by  de  Hoffmann 
and  Teller  [2]o   Here  we  only  mention  that  for  the  fast  perpendicular 
shock,  as  well  as  for  the  "non-Alfven"  parallel  shockp  the  jximp  of 
the  velocity  u  is  in  the  normal  direction^  as  may  be  seen  from 
(8o5)o''   Therefore,  if  in  these  cases  the  shock  is  observed  from 

The  slow  perpendicular  shock,  for  which  m  =  0,,  and  the  "Alfven" 
parallel  shock,  for  which  m  "   B  /[it  g    are   not  governed  by 
equations  ( 8ol4-)  =  (  806 ) ,  but  instead  by  equations  analogous  to 
{^ol\.)    and  (I4.06)  respectivelyo 
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an  appropriate  frame,  the  flow  velocity  will  be  norms.l  on  bouh 
sides, 

Fsr  "obliaue"  shocks.,  those  for  which  E  7^  B  n  .?^  0^  L.SLSz^ 
can  be  so  _ introduced  that  the  flow  veloci-!-.y  u  is  parallel  to  the 
magnetic  field  3  on  both  sides  of  the  front.   This  fact,  emphasized 
by  de  lioffmrmn  and  Teller,  can  be  read  off  from  relation  D^  ,   One 
need  only  choose  the  frame  such  that  u  =  B'  (u   -  U)B  on  one  side  J 
D..  then  implies  that  this  relation  also  holds  on  the  other  side. 
In  viev  of  (1)  the  electric  field  E  vanishes  in  this  frcrie  on  both, 
sides  cf  the  front. 

Using  formulas  {6.U)  to  (8,6,),  we  c&n  easily  describe  the 
.I'Timp  o:  the  abs_ol/ate_yalue  of__the  magnetic  field  across  tVxe  shock 
fron'.:.   ?rom  formula  IS,).!),,  we  have 

H e  n c ^ ,;,  u s  i  n  -7   f o i-mul  s  (8.6),, 

(6.7)      IB^]  ^   'i/[V]{  Tnr    '   r;J;/t.i)"^(3^  -  B^). 

Since  we  have  aBsum.sd  the  normal  n  to  point,  in  tho  direction 
in  which  the  fluid  cro.-'ses  the  shock  fronts  the  statoman:.  m-ide 
preceding  formula  (7.2'  iiriplies  that  [V]    <  0,   in  view  of  ( 8.3)  y 
we  are  therefore  able  to  state;   The  r^ap;netic  field  jLtillSriE.L'i  ^^: 
rist-s  a.crq ss  a  j'&st_shock  and  drops  across  _.a_slpv;  y-Jtock.   [■;•?- 
virtxie  of  D   the  pame  Is  true  of  the  magnituds-:  of  tb.e  tan^?ential 
component  of  B,   This  tnct  will  prove  T>articu"i, arly  s  igrj,ifi'.::;^.'nt  in 
connection  with  specific  flow  problems, 

Formula  (8.1^.)  implies  that  the  tangential  comc-onetit  of  the 
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magnetic  field  jumps  In  its  own  direction;  ioe„,  the  j\imp  [B] 

is  parallel  to  the  tangential  component  of  B   and  has  the  same 

o 

sense  as  the  tangential  component  of  Bo   Thus  across  a  fast  shock 


the  tangential  component  of  the  magnetic  field  retains  its  direction 
and  Increases  its  magnitude,  while  across  a  slow  shock  this  component 
may  retain  or  reverse  its  direction  and  must  decrease  in  magnitude o 

The  analysis  of  shocks  as  given  here  seems  appropriate  if  ©ne 
desires  to  obtain  a  quick  siorvey  of  the  possible  types  of  shocks 
and  to  derive  some  of  their  simple  propertieso   In  this  analysis 
we  have  extensively  used  mean  values,  involving  values  on  both 
sides  of  the  front o   In  a  niiraerical  problem,  however,  the  more 
important  questions  are  usually  those  which  refer  to  the  behavior 
of  the  various  quantities  on  each  side  separatelyo   To  answer  such 
questions,  the  analysis  of  Ericson  and  Bazer  [11]  may  be  m©re 
appropriate o   We  summarize  below  various  results  concerning  such 
problems,  the  proofs  of  which  may  be  found  in  [ll]o 

One  question  that  may  arise  is  whether  or  not  a  tangential 
component  of  the  magnetic  field  may  be  produced  through  a  shock 
if  it  was  absent  ahead  of  it,  or  whether  or  not  it  may  be  wiped 
out  if  it  was  present  ahead  of  it.   Shocks  through  which  this 
happens  may  be  called  "complete  switch-on"  or  "switch^eff "  shocks » 
Clearly,  a  complete  switch-on  shock  must  be  fast  or  a  complete 
switch~®ff  shock  must  be  slowo 

It  can  be  shown  that  complete  switch-on  shocks  exist  only 

if  the  Alfven  speed  is  supersonic  ahead  of  the  shock;  and  even 

then  only  if  the  shock  strength  (e<,go  measured  by  the  pressure 

ratio  Pt/p  )  lies  bel©w  a  critical  value o   If  one  lets  the  shock 
X   ®  — 
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strength  increase,  the  gain  in  tangential  component  of  B  first 
increases  and  then  decreases,  becoming  zero  when  the  critical 
stirength  is  reached.   On  the  other  hand,  complete  switch-off 
shocks  always  exist  if  the  Alfven  speed  behind  the  shock  front  is 
subsonic;   if  this  speed  is  supersonic  there,  they  exist  only  if 
the  shock  strength  exceeds  a  critical  value . 

Gas  dynamical  shocks  have  the  property  that  the  normal  component" 
of  the  flow  velocity  relative  to  the  shock  velocity,  u   ~  IK    is 
supersonic  ahead  of  the  shock  front,  i.e.  on  the  side  (0)  from  where 
the  fluid  comes,  and  subsonic  behind  it.   Por  a  hydromagnetic  shock 
we  may  state:   The  normal  flow  velocity  relative  to  a  fast  shock  is 
greater  than  the  fast  disturbance  speed  ahead  of  the  front  and  less 
than  it  but  greater  than  (or  equal  to)  the  transverse  distiirban.:e 
spee i  behind  iti 

u  -  U  >  c„   .      on  side  (0)^  ahead, 

b  <  u  -  U  <  c„   ,       on  side  (l),  behindo 
n  -        n        last 

The  e Duality  sign  holds  only  for  complete  switch-on  shocks.   The 
normal  flow  velocity  relative  to  a  slow  shock  is  greater  than  the 
slcv  dist^JTbance  speed  ahead  of  the  front ^  but  less  than  (or  equal 
to}  the  transverse  speed  on  both  sides  of  the  front ; 


u  -  U  <  b         on  side  (l),  behind, 
n     —  n  J: 

b   >  u   -  U  >  c  ,        on  side  (0),  ahead, 
n  -   n        slow  * 

I'he  second  equalit;/-  sign  holds  only  for  complete  switch-off  shocks. 
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These  facts  may  be  illustrated  In  Figures  2  and  2A,  in 
which  the  shock  is  assumed  to  be  stationary,  U  =  Oo   Only  those 
disturbance  motions  are  shown  that  travel  against  the  flowo 


Just  as  in  gas  dynamics,  a  fast  or  a  slow  hydromagnetlc  shock 

is  determined  by  prescribing  all  quantities  ahead  of  it  and  the 

pressiore  p  >  p    or  the  relative  velocity  v  =u   =U<v 
*- — — — '■—  — — ^  n    n        n 

behind  ito   It  is,  however,  not  possible  to  prescribe  arbitrarily 
the  tangential  component  of  the  magnetic  field  behind  the  shock 
front,  and  where  it  is  possible  to  prescribe  it  at  all,  the  shock 
may  not  be  iinlquely  determinedo   This  fact  is  clearly  indicated  by 
the  remarks  about  switch-on  and  switch-off  shocks  made  abovoo 

An  important  insight  into  the  connections  between  the  various 
types  of  shocks  may  be  obtained  by  consideration  of  the  family  of 
slow  shocks  starting  from  a  given  state  (with  a  non=zero  tangential 
component  of  B)  ahead  of  the  front o   As  the  shock  strength  increases, 
the  tangential  magnetic  field  component  behind  the  front  first 
decreases  to  zero,  then  Increases  with  direction  reversed,  approach- 
ing in  the  limit  the  negative  of  its  value  aheado   In  other  words, 
a  continuous  transition  may  be  made  from  a  weak  slow  shock,  such  as 
is  depicted  in  Figiire  2A,  through  a  complete  switch-off  shock  to  a 
special  case  of  the  "transverse"  shock  considered  in  the  following 
sectiono   This  remarkable  fact  was  first  noticed  by  Ericson  and 
Bazer  [11] o 
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9o   Transverse  Shocks,   Contact  Discontinuities 
The  root 

(9.1)  m^^   =  (B2/^.r)^/2  ^ 

of  equation  det  (E)  =  0,  cf,  (7o7)j  is  associated  with  a  type  of 
shcrCk  which  will  be  called  transverseo"   Such  transverse  shocks 
correspond  to  the  transverse  waves  discussed  in  Section  I;.;  they 
are  sometimes  called  "Alfv^n  shocks."   In  analogy  to  expressions 
(l4.,2)  we  find  the  expressions 

(9.2)  [B]  =  kmB  /  n, 

(9.3)  [u]  =  l?;i"-'-B^B  X  n, 

(9.I4-)      [t^]    =  0, 

(9.5)  [S]  =  0. 

The  last  relation  is  derived  from  (7.2)  and  [t]   =  0.   Also, 
[u  ]  =  0  holds,  as  seen  from  E^.   Thus,  the  only  quantities  that 
.•jump  across  a  transverse  shock  are  the  tangential  components  of 
the  magnetic  field  and  of  the  velocity. 

Relation  (9.2)  implies  relation 

(9.6)  [3^3   =  0. 

In  other  words,  the  strength  of  the  magnetic  field  is  unchanged 
across  a  transverse  shock.   The  magnetic  field  therefore  rotates 
in  the  plane  of  the  shock  and  the  flow  velocity  undergoes  a 
tangential  change  parallel  to  [B],   All  other  quantities  remain 
continuous.   The  possibility  of  this  occtirrence  necessitated  the     J 

qualification  made  in  the  statement  preceding  formula  (7.2)0 

—  -— _— ^  ^  ___         - 

A  transverse  discontinuity  is  called  a  "shock"  since  fluid  crosses 
the  STorfacei  in  other  respects  it  is  closer  to  a  contact  discontinue tj 
See  [16], 
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In  general,  the  angle  through  which  the  tangential  magnetic 
field  B  =  B  -  B  n  rotates  in  crossing  a  transverse  shock  is 
arbitrarjo   An  interesting  special  case  is  that  for  which  the 
rotation  angle  is  l80  ,    so  that  the  tangential  magnetic  field 
behind  the  shock  is  directed  oppositely  to  that  in  front o   As 
pointed  out   in  Section  8,  this  special  case  also  occurs  as  the 
limit  of  a  slow  shock  when  the  shock  strength  approaches  its 
maximum  admissible  value »   Thus  transverse  shocks  are  continuously 
connected  with  the  slow  shocks  considered  earlier,.   This  connection 
is  more  clearly  shown  in  Figxrre  3o 

Ret\irning  to  the  general  transverse  shock,  it  is  obvious  from 
(9o3)  that  a  frame  can  be  introduced  such  that  the  flow  is  parallel 
to  the  front  on  both  sides o   Therefore,  the  term  "transverse" 
shockj  introduced  by  de  Hoffmann  and  Teller j,  would  seem  appropriate o 

The  possibility  of  such  a  transverse  shock  is^  of  course^  to 
be  understood  as  a  mathematical  possibility^  referring  to  the 
existence  of  solutions  of  the  equations  (E)  for  the  root  (9ol)o 
Whether  or  not  such  shocks  are  possible  in  nature  is  another 
questioHo   In  factjj  it  would  seem  that  they  are  possible  only 
vmder  vmusual  circiimstances|  cfo  Section  12, 

The  root  m  =  0  of  equation  (7o7)  corresponds  to  a  discontinuity 
of  the  normal  component  B  of  B.   However,  this  possibility  is 
excluded  by  condition  D  5  cfo  the  arguments  in  Section  5»  which 
could  be  carried  over  herco 

A  second  root  m  =  0  would  have  occurred  if  we  had  not  omitted 
relation  D,  in  changing  the  system  (D)  over  into  the  system  (E)o 
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The  corresponding  discontinuity  would  be  a  proper  contact  discon- 
tinuity involving  no  flow  across  the  front,  i.e. 

u^  -  U  =  0. 

In  case  B  ^0  relations  (D)  then  imply 

[u]  =  0,   [B]  =  0,   [p]  =  0. 

The  latter  relation  does  not  require  [p]  =  0,  since  [S]  ^  0  is 
compatible  with  Di  .   Because  of  [B]  =  0  we  have  a  piirely  gas 
dynamical  contact  discontinuity,  in  fact,  since  [u]  =  0,  a 
special  one. 

In  case  B   =0,  on  the  other  hand,  we  can  only  conclude 

[p  +  bV2m-]  =  0, 

while  the  tangential  components  of  u  and  of  B  may  undergo  any 

jiomps.   In  fact,  the  contact  discontinuities  corresponding  to 

both  roots  ra  =  0  coalesce  in  the  case  B  =  0, 

n 

A  contact  discontinuity  which  involves  a  jump  in  the  tangential 
flow  velocity  will  be  called  a  "shear  flow  discontinuity."   It  is 
remarkable  that  in  a  conducting  fluid  no  shear  flow  discontinuity 
can  be  maintained  if  the  magnetic  field  has  a  normal  component 
B  f^  0  at  the  front.   This  fact  will  be  the  starting  point  for 
our  discussion  of  special  f].ow  problems  in  Section  13, 
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10 o   Simple  Waves" 


Naturallys,  the  hydro  =magne tic  flows  most  easily  accessible 
to  treatment  are  the  one=>dimensional  flows »   These  are  characterized 
by  the  condition  that  all  quantities  depend  only  on  one  space  variable; 
X  say  J,  in  addition  to  the  time,  and  hence  are  constant  on  each 
(yp2)°plane  at  each  time.   No  restriction  need  be  imposed  as  to 
the  presence  or  absence  of  the  y  and  z-components  of  the  vectors 
B  and  Uo 

The  problem  of  one-dimensional  waves  in  a  compressible  con- 
ducting medium  is  certainly  not  the  most  urgent  problem  of  magneto-= 
hydrod3mamics  that  needs  to  be  solved|  nevertheless ^  as  in  gas 
dynamics  J  the  study  of  such  problems  contributes  to  an  iinderstand- 
Ing  of  significant  hydro-magnetic  phenomenao 

In  gas  dynamics^,  the  simplest  types  of  one=dimensional  flows 
are  the  30°called  "simple  waves" o   Because  a  flow  region  adjacent 
to  a  region  of  constant  state  is  always  a  simple  wave,  these  waves 
may  be  used  very  effectively  as  building  blocks  in  constructing 
solutions  of  flow  problems,,   (A  state  of  flow  is  referred  to  as 
constant  in  a  region  if  all  significant  flow  quantities  are  time- 
and  space-independent  in  the  regiono ) 

Simple  waves  are  also  possible  in  magneto-hydrodynamicso 
They  have  essentially  the  same  properties  as  those  in  gas  dynamicso 
Moreover;,  it  appears  that  they  also  can  be  effectively  used  as 
building  blocks o 

A  simple  wave  may  be  characterized  by  describingj  not  the 

For  a  theory  of  simple  waves  associated  with  general  systems  of 
differential  equations  whose  coefficients  do  not  involve  the 
independent  variables,  see  Po  Lax  [l6]o 
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motion  of  a  particle,  but  the  motion  of  a  "phase" ^  as  determined 
by  a  set  of  values  for  the  quantities  B,  u^  po   In  a  simple  wave 
each  phase  moves  with  constant  velocity„   That  is,  the  values  of  B, 
Uj  p  are  constant  on  certain  straight  lines 

(lOol)     X   =  Ut  +   C 

in  the  {x,t)  plane »  The  number  U,  as  vjell  as  B,  u,  and  p,  may  be 
considered  a  fiinction  of  the  parameter  £,«  Instead  of  U  we  intro- 
duce the  quantity  c  >  0  defined  by 

(10„2)     U  =  u^  +  0. 

This  quantity  c  will  turn  out  to  be  one  of  the  characteristic 
distiirbance  speeds?  i.e.,  one  of  the  roots  of  equation  {2o6), 

From  this  characterization  of  simple  waves  one  may  derive 
the  following  recipe  for  setting  up  the  equations  governing  them^ 
In  equations  A  (see  Section  1),  one  should  replace  V  and  -^'r   by 
■^d/d^s,  0,  Of  and  -{u  t   cjd/d^p  oPp  leaving  open  the  choice  of 
the  parameter  describing  the  phase,,  simply  by  ^  dp  0,  0 '^  and 

=  (ti^  ±  C)do 

The  equations  thus  obtained  are 


F^        dB^  =  0; 
o  X      ' 


P,        xc  dB  +  B  du   =>  B  du  =  0, 
1        +    y    y  X    X  y     * 


TC  dB„  +  B„du  =  B^du  -  0, 
+    2     z   X    X   z     ' 
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P„        IPcdu  +  dp  +  d(B  ,/2n) 


y   "   X  y 


-pcdu^  +  \i.     B^dB^  =-  0, 


■=1 
+pcdUg  -  ti  B^dBg  =  0, 

P3        ^cdp   +  pdu^  =  0, 

P^        ;cdS  =  0„ 

When  the  equations  F-  to  P.  are  considered  linear  equations 
for  the  seven  differentials  dB  , , o o ,  dS   it  is  clear  that  the 
determinant  of  this  system  must  vanisho   One  verifies  ■=  by  direct 
computation  or  by  comparison  with  equations  (B)  in  Section  2  =»  that, 
except  for  a  single  factor  of  ^c,  this  determinant  is  precisely 
det  (B),  cfo  (2o5)s  with  the  normal  direction  n  taken  as  the 
x-directioHo   Hence,  the  speed  c  must  be  one  of  the  roots  of  the 
equation  det  (B)  -  Oj,  ioeo?  one  of  the  characteristic  distxarbance 
speeds  corresponding  to  the  phase  B,  u^  Po 

According  to  which  kind  of  speed  c  enters,  the  simple  wave 
is  fastg  slow J  or  transverse |  if  c  =  0  the  wave  is  a  "contact 
layer" o 

As  we  shall  see  in  the  subsequent  sectlonsj,  the  differential 
equations  (F)  can  be  greatly  simplifiedo   In  the  case  of  a  poly° 
tropic  gaSp  they  can  even  be  solved  explicitlyo 
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llo   Fast  and  Slow  Simple  Waves 

Equation  (3<.2)  for  the  fast  and  slow  disturbance  speeds  may 
in  the  present  case  be  written  in  the  form 

(11.1)  ti""^(B^  +  bI)    =  c'^Cc^  -   a.^){pc^    -  B^/tt), 

y     z  jt 

where,  cfo  (2), 

a  =  dp/dp o 

It  is  more  convenient  for  the  present  purpose  to  introduce  the 
square  ratio 

(11.2)  q  =  c^a^ 

of  disturbance  speed  to  sound  speed  as  independent  variable  and 
to  try  to  express  all  other  quantities  in  terms  of  q»   It  is  also 
convenient  to  introduce  the  square  ratio 

(llo  3)      s  =  a^b^  =  ^ipa^/B^ 

of  sound  speed  to  "normal"  Alfven  speed  as  dependent  variable,, 
Relation  (11,1)  then  becomes 

(11, U)     B^  +  B^  =  (q  -  l)(s  -  q°^)B^„ 

Prom  relations  P^        and  Po   one   fiirther  derives   the  relation 

(11.5)  2(q    -    Dds    =  rd(q    =    DCs    -    q"^), 

where 

(11.6)  r  =  l+P^(log||). 

(See  the  Appendix  for  a  more  systematic  derivation  of  this  relation.) 
For  polytropic  gases^  y  is  a  constant,  and  relation  (11,5) 
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becomes  a  linear  differential  equation  for  s  as  a  function  of  q, 
which  can  be  solved  explicitly,,   Prom  relations  F.  and  Pp   ,  Pp 
one  infers  that  the  ratio  B  /b  is  constant^  so  thatj,  except  for 
this  constant^  B   and  B   can  be  determined  explicitly  in  terms  of 
the  constant  B  by  (lloU.)o   Since  p  and  a  may  be  regarded  as  knovm 
functions  of  s^  cf„  (llo3),  also  c  =  a/q  may  be  considered  knowno 
The  velocity  u  can  now  be  determined  by  integrating  the  differentials 

(llo7)      du^  =  tcp°^dp, 

(lie.  8)      dUy  =  ;(npc)"^B^dBy, 

du^  -  :f(tipc)°^B^dB^„ 

As   implied  by  relation   {3ol\.)g 

2  2  =1 

(llo9)  |ipc     >  B        or   s   >  q  in  a  fast  wave^, 

2    2  ~1 

p,pc  <  B   or  s  <  q    in  a  slow  wave^ 

thus  the  wave  speed  c  may  agree  with  the  Alfven  speedo   On  the 
other  hand,  it  could  be  shown  that  the  wave  speed  c  can  never 
coalesce  with  the  sound  speed  aj  ioe,, 

(llolO)     q  >  1   in  a  fast  wave^ 

q  <  1   in  a  slow  wavoo 

In  a  fast  wave  the  compression  may  tend  to  infinity,  and  hence 
s  °>  oo,  while  in  a  slow  wave  cavitation  may  be  reached^,  s  =>  Oo 
Of  course  J,  a  fast  (or  slow)  wave  may  be  either  a  compression 
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wave  or  a  rarefaction  waveo   The  tangential  components  of  the 
magnetic  field  |B  [  and  JB  i  increase  across  a  fast^  decrease 
across  a  slow  compression  wave,  but  they  decrease  across  a  fast,, 
and  increase  across  a  slow  rarefaction  waveo 

Finally,  we  mention  that  "centerei  rarefaction  waves"  exist, 
as  in  gas  dynamics.   They  are  characterized  by  the  condition  that 
at  some  initial  time,  say  t  =  0,  all  phases  involved  are  concen- 
trated at  the  same  point.   (One  then  must  set  £,  =   const,,  in  formula 
(10„1)  and  use  a  different  parameter  to  characterize  phases,)   If 
we  imagine  that  a  centered  wave  separates  tv;o  constant  states  of 
the  fluid,  these  states  will  be  adjacent  at  the  time  t  =  0,   In 
other  words,  the  centered  wave  then  resolves  an  initial  discon- 
tinuity , 
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12 o   Transverse  Waves  and  Contact  La'Vft'^s 

While  the  ratio  of  H^  to  H   is  constant  across  fast  and 
slow  waves y  this  is  not  the  case  for  the  transverse  simple  vaves^ 
which  we  are  now  going  to  describe.   The  speed  c  of  these  waves 
is  the  Alfven  speed 

(12.1)    c  =  (b;/^p)^/2, 

Piirthermore^  B  .  p,  and  S,  and  therefore  p  and  ii  .  are  constant 
across  these  waves.   The  tangential  magnetic  field,  hoivever, 
rotates: 

(12o2)      B  =  G  cos  0,   B^^,  =  u  sin  ©, 

with  0  -  0(x}  being  any  function  of  the  phase »   B'lu^ther, 

(12.3)       u  ==  Q.  +  gB'-^B  ,   u  =  a   7  cB-'^B  , 

with  any  numbers  a^,  a  „   Thus^  the  flow  in  a  transverse  wave  may 
be  considered  a  shear  flow.   All  particles  on  the  same  {y,2)-plane 
move  in  the  same  straight  line.   This  shear  flow  is  evidently  a 
steady  flow  if  it  is  observed  from  a  frame  with  respect  to  which 
u  1  c  =  0,   Thus,  in  contrast  to  the  situation  in  gas  dynamics, 

•A* 

there  do  exist  non-constant  stea dy  flows  in  magneto-hydrod-Tnamlcs „ 
In  order  to  maintain  such  a  shear  flow,  it  would  be  necessary 
to  have  at  large  (y, z)~distances  a  mechanism  which  supplies  the 
velocity  u  in  the  proper  directions  there.   It  would  seem  that 
such  a  mechanism  would  not  easily  occur  under  natural  circumstances 
and  that  it  would  have  to  be  rather  artificial. 
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A  transverse  wave  may  connect  two  constant  states  vjith 
different  velocities  u  and  magnetic  fields  B,   One  may  Imagine 
the  layer  covered  by  the  wave  to  be  arbitrarily  thing  one  then 
may  approximate  the  transition  by  a  discontinuity.   This  discon- 
tinuity would  be  exactly  a  transverse  shock.   ThuSe,  a_  transverse 
shock  may  be  considered  the  limit  of  transverse  simple  waveso   The 
remark  about  the  artificial  nature  of  transverse  waves,  therefore,, 
apply  just  as  well  to  transverse  shockso 

Only  a  short  remark  need  be  made  about  "contact  layers".   In 
accordance  with  c  =  0^  there  is  no  flow  across  such  a  layers:  1,6 o;, 
u  is  constant  across  it.   If  B  ~   Oj,  all  other  q\iantities  may 
vary  except  that  the  relation 

(12.il.)       p  +  (B^  +  B^)/2|i  -  constant: 

should  hold.   If  B  ^  0,,  however,  B,  u,  and  p  are  consfcanf; .,  and 
only  p  and  S  may  vary»   Thus,  in  contrast  to  gas  dynamics ,  a  shear 

flow  layer  across  which  the  tangential  flow  components  u   ^    u  .  vary 

y   2 

cannot  be  maintained  in  a  conducting  fluid  if  the  magnetic  field 
possesses  a  normal  component  B  ^  0, 
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13.   The  Resolution  of  a  She^r  Plcw  DiscontinuitT 

A-:-  was  explained  at  the  end  of  Section  9,  no  shetr:--  flow  dis- 
continuity can  remain  unchanged  in  tne  presence.:  of  a  no.gnetic  field 
whose  rorinal  coriDonent  is  not  ^ero.   If  ?!uch  a  discontinuity  is 
present  s.t   an   initial  time,  a  wave  motion  must  penult  xvhich  resolves 
it,   Thf:  nature  of  this  wave  motion  todll  now  be  described. 

Specifically,  we  consider  the  following  problem,   /..t  an  initial 
time,  t  =  C,  the  fluid  is  at  rest  on  one  side,  x  >   0,  of  the  plane 
X  =  0,  1,'hlle  on  the  ether  side,  x  <  0_,  it  possesses  a  constant 
tangential  velocity  (u,..,  u  )  ^   0^  but  no  normal  veloci-cy,  so  that 
u^  =  0,   The  density  is  constant  and  the  sar.e  on  both  sides,  and  a 
constant  magnetic  field  is  present  -  the  same  on  both  sides  -  with 
a  non- '/a '.lishing  normal  component  B  ^  0. 

The  impossibility  of  maintaining  a  shear  flovr  discontinuity 
may  be  visualized  as  follows.   Instead  of  a  discontinuity,  consider 
a  thin  shear  floxv  layer  across  which  a  tangential  flow  component, 
say  u  .  varies  smoothly  from  a  positive  value  on  the  left  hand  side 
to  the  value  zero  on  the  right;  for  slm.pl lei ty  assume  u^  =  0 
throughout.   Also  assume  E  ~-  3^  =  0  and  B^  >  0.   ^'he  basic  assump- 
tion of  magneto-m'-drodynamlcs,  em.bodled  in  formula  (1),  now  implies 
that  the  electric  field  component  S   does  not  vanish  on  the  left 
but  does  so  on  the  right.   Therefore,  curl  E  =  0  in  the  layer. 
From  A-,  \-ie    then  conclude  that  the  tangential  magnetic  field, 
specifically  the  component  -B  ,  will  grow  within  the  lajrer.   Since 
this  component  at  first  remains  zero  outside  of  the  layer,  its 
curl,  and  hence  the  current  J,  will  be  different  from  zero  in  the 
layer.   Specifically,  the  component  J  wiT  1  vi-ry  from  negative  to 
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positive  values  across  the  layer«  Hence  the  fluid  particles  in 
the  layer  will  experience  a  force  whose  x-component  varies  from 
negative  to  positive  values  across  the  layer  and  hence  tends  to 
push  the  particles  away  from  the  layero 

Before  describing  the  details  of  the  resulting  wave  motion  we 
mention  a  concrete  situation  in  which  such  initial  shear  flow 
discontinuity  may  occ\ir|  namely,  if  a  jet  of  conducting  fluid 
shoots  into  conducting  fluid  at  rest  so  quickly  that  the  hydro- 
magnetic  adjustments  just  discussed  have  not  yet  developed,   A 
similar  case  was  described  by  Alfven  [1]  in  explaining  the  possible 
origin  of  hydro-magnetic  waves.   Assuming  the  jet  to  proceed  in  the 
y-direction  and  to  be  much  wider  In  the  z-direction  than  in  the 
x-direction,  we  may  approximate  it  by  a  constant  flow  in  the  y- 
direction  between  two  parallel  planes  x  =  constant.   The  waves  which 
result  from  the  two  interfaces  will  interact  only  after  some  time. 
Up  to  this  time,  therefore,  the  situation  may  be  described  in  terms 
of  the  wave  motion  resulting  from  a  single  interface » 

We  maintain  that  the  resolution  of  the  shear  flow  discontinuity 
is  effected  by  two  fast  shocks  followed  by  two  slow  centered  rare-' 
faction  wavesc   After  the  waves  have  formedj,  the  fluid  has  acquired 
the  mean  tangential  velocity,  while  a  tangential  component  of  the 
magnetic  field  has  increased  (or  developed  if  none  was  present 
originally), 

A  diagram  of  the  resolving  flowp  in  which  only  the  x-compon- 
ent of  the  particle  motion  is  indicated,  is  given  in  Figure  I4., 

It  may  be  mentioned,  incidentally,  without  giving  supporting 
arg\aments,  that  after  the  waves  coming  from  the  two  interfaces 
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have  interacted  the  fluid  will  come  to  rest  in  the  jet  region 
while  the  tangential  magnetic  field  has  a  further  increased  value, 
at  least  for  some  timeo 

In  attacking  the  problem  in  detail^  we  have,  for  simplicity^ 
assToraed  that  all  z-components  vanish  and  that^  at  time  t  =  Op  B 
is  a  positive  constant  and  B  =  Oo   Fxirther,  we  have  assumed  that 

y 

p  and  p  are  constant,  the  saine  on  both  sides, 
(l3ol)       P  =  P,   P  =  P, 

and  that  u  -  0  on  both  sldeSo   For  convenience  we  have  assumed 
that  the  flow  is  observed  from  a  frame  moving  with  the  mean 
tangential  velocity^,  so  that 


X  >  0, 


(13o2)      u 

J  I  =11   ■»"  <j?*  0 


o 

u  being  a  positive  constant, 


The  conditions  to  be  satisfied  after  the  passage  of  the 
waves  (ioeo,  in  the  region  of  constant  state  containing  x  =  0) 
simplify  because  of  the  symmetry  of  the  problemo   They  are 


(13.3)^     u^     ^     0 


s 


(13.3)      u   =  0, 

Since  B   is  even  in  x„  no  condition  need  be  imposed  on  B  o 
7  '  1 

It  was  mentioned  at  the  end  of  Section  8  that  a  shock  would 
be  determined  if  one  quantity  such  as  p  or  u  is  prescribed  behind 
it,  provided  the  state  in  front  of  it  is  known.   Except  for  a 
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qualification  to  be  discussed  below,  the  same  is  true  for  a  simple 
waveo   It  is  J,  therefore,  natural  to  expect  that  two  quantities  could 
be  prescribed  behind  a  pair  of  waves  when  the  state  in  front  is 
known.   It  is  not  obvious,  however,  whether  or  not  the  two  velocity 
components  u  ,  u  may  be  prescribed  without  limitation  behind  a 
pair  of  waves.   In  the  present  problem  this  appears  to  be  the  case. 

The  qualification  mentioned  above  is  this;   Suppose  a  piston 
at  one  end  of  a  gas  filled  tube  is  withdrawn  with  speed  greater 
than  a  certain  "escape  speed",  cf,  [1^],   Then  the  resulting 
rarefaction  wave  will  lead  to  "cavitation".   The  piston  will 
separate  from  the  gas  and  a  vacuxim  zone  will  be  formed.   The  gas 
at  the  edge  of  this  zone  will  move  with  the  escape  velocity  and 
not  with  the  piston  velocity,   A  similar  phenomenon  may  occior  here, 
so  that  condition  (13. 3)^  must  be  modified.   It  should  read 

(I3,i|)^     Either  u^  =  0  or   p  =  0, 

thtis  permitting  the  presence  behind^the  rarefaction  waves  of  a 
vacuum  zone  which  expands  with  an  appropriate  escape  speed.   When 
this  vacuum  zone  is  present j.  condition  (13,3)   must  be  replaced 
by  the  condition  that 

(I3»l4.)       u  B   -  u  B   =0 
^  ^'y      y  X    X  y 

at  the  edges  of  the  vacuum  zone.   As  noticed  by  Bazer  [10],  this 
is  the  correct  condition  to  insiire  that  the  electric  field  E  =  B/u 
vanishes  there  and  In  the  whole  vacuum  region,  ^i. 

For  the  description  of  the  wave  motion  in  detail  it  is 
necessary  to  solve  a  number  of  transcendental  equations  which  are 
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expressed  in  terms  of  explicitly  given  integrals  derived  from 
the  considerations  of  Section  11» 

Of  particular  interest  is  the  magnitude  of  the  tangential 
component  B   of  the  magnetic  field  that  has  developed  at  the 
center  x  =  0  after  passage  of  the  waves.   For  small  values  of 
the  initial  shear  flow  discontinuity  2u  ,  this  component  may  be 
described  by  the  formula 

(13.5)        B      =  (TCtip)^/^  S     , 

o 

in  which    X    in  t\irn  depends   on  u   . 

y 

If  u  approaches  zero,  X.  approaches  the  value  lo   In  fact, 
this  same  value  of  X.    would  have  been  obtained  if  one  had  assumed 
the  fluid  to  be  incompressible  and  had  described  the  resolution 
of  the  discontinuity  with  the  aid  of  two  Alfven  waves,  one  traveling 

in  each  direction, 

o 
If  u   is  sufficiently  large,  cavitation  occurs.   In  such  a  case, 

the  soTind  speed  eventually  becomes  less  than  the  Alfven  speed  even 

if  it  was  much  larger  originally.   Thus  the  medium  could  certainly 

not  be  considered  incompressible.   The  limiting  form  of  the 

o 
tangential  component  B   as  u  ~>  oo  is  given  by 

(i3o6)    By  -  [U+r)\iP^]^^^  (B^Sy)^/2^ 

where  y  Is  as  defined  by  (11,6). 

The  derivation  of  formulas  (13.5)  and  (13«6),  as  well  as 
complete  details  of  the  resolving  flow,  are  given  by  Bazer  [lOl, 
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APPENDIX:   Simple  Waves  and  Groups  of  Trans  format,  ions 
by  K,  von  Hagenow 

In  the  following  It  is  demonstrated  how  the  IrJierent  sytfimetry 
of  the  equations  of  one-dimensional  flow  can  be  ufei  to  achieve 
the  reduction  of  the  equations  of  hydromagnetic  simple  waves  (see 
Section  10)  to  a  siimple  ordinary  differential  equation. 

Simple  waves  are  special  solutions  of  the  eauations  of  one- 
dimensional  flow  which  are  obtained  from  equations  A^  through  A, 

1       ^        k 

of  Section  1  by  assuming  all  quantities  to  depend  only  on  the 
coordinate  x  and  the  time  t.   The  divergence  condition  /■«  then 
reduces  to  B^  =  const.   The  equations  are  non-rclatl vistie ,  End, 
as  physically  no  reference  S7/-stem  can  be  preferred,  they  must  be 
invariant  under  the  ^-parameter  group  consisting  of  the  following 
transformations:   a  translation  with  constant  speed  in  an  arbitrary 
direction  and  a  rotation  arotmd  the  x-axis  by  an  arbitrary  anp-le. 
The  direction  of  the  x-axis  has  to  remain  unchanged,  of  course,  for 
we  have  distinguished  it  from  the  other  two  coordinate  axes  by  the 
assiimption  of  one-dlnensional  flow,    Polytropic  gases,,  where  the 
pressure  p  is  given  by 

(ll+.l)  p  =  p^'  e^ 

allow  the  additional  similarity  transformation 


X   =  ax  , 


(1I|.,2}         t'  =  t  . 


vjith 
(l)+.,3) 


u'     =  a  u 


S    =^  S  +  2y  log  a 
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and  pressure  p  and  magnetic  field  B  unchangf-d.   If  D  =0,  the 
pressure  p  and  field  B  can  also  be  scaled  by  a  factor: 

P'   =  Pp 

Hk^h)  ,     -,/p 

B   =  p  /  B. 

Now  in  simple  t>faves,  all  variables  u,B,p,S  are  constant  on  a 
family  of  straight  characteristics   (of.  eq.  10,1) ,   As  our  trans- 
formations map  straight  lines  into  straight  lines  (in  the  x-t  plane) 
and  as  invariance  of  the  eouation  implies  mapping  of  characteristics 
into  characteristics,  they  also  map  simple  \\raves  into  simple  waves. 
That  is,  the  system  of  ordinary  differential  equations   P-  through 
P,  (section  10)  is  invariant  under  the  corresponding  transformation 
of  the  dependent  variables. 

Let  us  introduce  new  variables  as  follows:   the  pressure  p, 
the  square  of  the  field,  B  ,  the  angle  0   between  the  y-z  field 
component  and  the  y-axis,  the  velocity  components  u  in  the  x- 
direction,  V  in  the  direction  of  above  field-component,  W 
orthogonal  to  V  in  the  y-z  plane,  and  the  entropy  S.   Then,  if  we 
write  the  system  P-,  through  Pi  of  Section  10  symbolically: 

(11^.5)  ^  a^^   dwj^  =(y^, 

we  see  that  we  can  add  an  arbitrary  constant  to  each  of  the 
variables  0,  u  ,  V,  W,  i.e.,  except  for  a  factor  common  to  each 
rovj,  the  matrix  a    can  but  depend  on  the  remaining  variables 
p,  B^,  S. , 
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S  is  constant  according  to  F,    so  the  system  tnust  be  reducible 
to  a  single  equation  involving  the  quantities  p  and  B  onlyj 
containing  S  as  a  parameter. 

Now  that  we  know  vjhat  equation  to  expect,  we  can  obtain 
the  reduction  without  actually  introducing  our  new  variables. 
Introducing  the  sound  speed  a  with 

(II4-.6)  dp  =   a  dp, 

and  using  equation  P-.,  we  can  write  equation  P^ 

(11^.7)  (l-q)dp  +  d(B^)/2ti   =   0 

with 

op 

(III.. 8)  q  =  c^/a         (cf,  equation  11.2) 


2 

end  we  know  that  q  depends  on  B  ,  p  and  S  only,  the  latter  being 

2 

constant  in  a  simple  wave  anyway.   We  can  therefore  express  B 

a  function  of  p  and  q,  and  ?et  exactly  equation  (11,5),  for  the 
quantity  s  defined  in  e^'uation  (11,3)  is  related  to  p  by 


as 


(ll+,9)  ydp  =  B^VM-ds 

with  Y  defined  by  equation  (11,6),  while  equation  (ll.h.)  allows 

2 
one  to  express  B   in  terns  of  s  and  q  for  fast  and  slow  waves. 

2 

Transverse  waves  are  trivial,  because  In  them  p  end  therefore  B 

are  constants,   Polytropic  gases  allovj  f-urthermore  the  transforma- 
tion (1)4,3),  l,e.  the  entropy  S  does  not  appear  in  (lk,7).   The 
special  case  of  vanishing  B   allows  the  transformation  (l[|,l4-). 
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2 

i.e.  q  mu3t  be  homogeneous  of  degree  o  in  p  and  B  .   It  is  note- 

vjorthy  that  we  got  our  results  without  making  explicit  use  of 

(11. I4.),  which  deter.Tiines  the  distiirbance  speeds.   If  we  actually 

2 
compute  B   in  terms  of  p  and  q  (for  polytropic  gases)  we    find  that 

it  is  linear  in  p,  i.e.,  equation  (llj..?)  can  be  integrated 

explicitly.   This  feet  cannot  be  deduced  by  looking  at  the 

symmetries  of  the  one-dimensional  equations,  but  we  can  make  the 

following  statement.   If  we  treat  also  ^x  as  sn  unknown,  then  the 

eauations  A   through  A.  allow  (ll|,2, 3,i^-)  even  for  B   7^  0 .   But  now, 

Dior  reduction  to  (Ik, 7)  allows  us  only  to  conclude  that 

dU.lO)         B^  =  B^(p,q,B^). 

The  vjave  formation  leaves  q  fixed,  for,  according  to  (l[i,2),  the 
characteristic  direction,  i.e,  c,  is  multiplied  by  a,  and  from 
(ll].,l)  we  deduce  that  the  snimd  speed  a  is  transformed  by  the  same 
factor.   Then  {l}x,L\.]    shows  that  B^  must  be  homogeneous  of  first 
degree  in  p  and  B   but  linearity  does  not  follow.   If  B  7^  0  it 
is  only  the  simplest  possibility  compatible  vjith  the  transformation 
properties  of  the  eauations,   Eauation  (Ik,?)  allows  an  Interesting 
interpretation:   With  (ll4.,6,8)  we  can  write  it 

2 
(1I|,11)       d(p  +  I-  )  =  c^dp 

i,e.  the  variation  of  the  total  "press\are""  with  the  density  gives 
the  square  of  the  disturbance  speed.   This  is  of  course  true  for  any 
longitudinal  vjave  in  a  compressible  medium. 


"We  should  really  talk  about  a  oressvire-tensor,  see  L\!ist  [?],  for 
the  magnetic  force  is  enisotropic,  but  for  one-dimensional  flow 
only  the  1-1  component  gives  a  force  in  the  x-directlon. 
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Thus  we  can  already  conclude  from  this  remark  that  the  fast  and 
slow  -        3  purely  longitudinal,  because  the  speed  of  transverse 
waves  is  in  general  not  given  by  (ll4.,ll),  but  depends  on  restoring 
forces  against  torsiono 
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FIGURE  1 
The  Three  Types  of  Disttirbance  Waves. 


Transverse 


B 


FIGURE  lA 

Wave  fronts  at  a  time  t  >  0  originating  from  a  point 
disturbance  at  the  time  t  =  0.  Their  envelopes,  also 
shown,  are  given  by  the  intersection  of  the  character- 
istic cone  with  t  =  constant  and  z  =  constant  (B„=  0), 
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FIGURE  2 


Fast 


Transverse 


Slow 


FIGURE  2A 
Transitions  Through  Stationary  Shock  Fronts 
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Transverse 


Transverse 


FIGURE  3 

Locus  of  the  possible  tangential  components 
of  magnetic  field  behind  the  various  types  of 
shock  corresponding  to  the  fixed  tangential 
component  B  ahead  of  the  shock. 
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Particle  Path 


Slow  Rarefaction  Wave 


Phase  Motion 


Past  Shock 


Past  Shock 
(Pure  gas -dynamical 
or  switch-on) 


Resolution  of  a  Shear  Plow  Discontinuity 


P 


-►  X 


Typical  Distribution  of  Pressure  at  Time  t,  >  0 


B 


Switch-on  shock 


■^  X 


gas  shock 


Typical  Distribution  of  Tangential  Magnetic  Field  at  Time  t^^  >  0 

FIGURE  k 
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